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The theoretical statements about a restoration of a superconductivity at magnetic 
fields higher than the quasiclassical upper critical field and a reentrance of superconduc- 
tivity at temperatures T C (H) » T c (0) in the superconductors with open Fermi surfaces 
are reinvestigated taking into account a scattering of quasiparticles on the impurities. 

The system of integral equations for determination of the upper critical field parallel to 
the conducting planes in a layered conventional and unconventional superconductors with 
impurities are derived. The H C 2(T) values for the "clean" case in the Ginzburg-Landau 
regime and at any temperature in the "dirty" case are found analytically. The upper limit 
of the superconductor purity when the upper critical field definately has a finite value is 
established. 
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I. INTRODUCTION 



The theory of upper critical field in highly anisotropic quasi-two-dimensional superconductors for the 
field orientation parallel to conducting layers, has been developed by A.Lebed' and K.YamajiEI. ft was 
shown that like in quasi-one-dimentional caseo at low enough temperatures (see below) the upper critical 
field starts to diverge such that a superconductivity is conserved in arbitrary large magnetic fields. 
Moreover in high enough magnetic fields (see below) the critical temperature of a superconducting phase 
transition has a tendency to restore of its zero field value. These statements are literally true for the 
superconductors with triplet pairing. For a singlet pairing the existance of a superconductivity in high 
fields is restricted by the paramagnetic limit H p . So, the low temperature stability of a superconducting 
state under magnetic field can serve as the indication on the superconductivity with triplet pairing. 

The tendency for low temperature divergence of the upper critical field H C 2 > H p in ..nuasi-two- 
dimentional K-type ET organic superconductors has been reported recently by T.Ishigurcd. Similar 
behavior has been observed earlier in quasi-one-dimensional organic superconductor (T MT S F)2R£f by 
the group of P.ChaikinQ. So, the observations being in correspondence with theoretical predictionsliro say 
in favor of triplet type of superconductivity in both type of these materials. 

Another popular layered superconductor is Sr2RuC>4. It demonstrates the properties compatible with 
triplet superconductivityD. Unlike isostructural high-T c cuprates La2- x Ba x CuO^ the normal state of 
layered perovskite oxide S^RuO^ conforms with the predictions of Fermi- liquid theoryu. Unconventional 
nature of the superconducting state in this material manifests itself through the strong supression of T c 
by nonmagnetic impurities&El as well as by the exibition of a sharp decrease without _the coherence peak 
of 101 Ru nucleare spin-relaxation rate 1/Ti followed by T 3 behavior down to 0.15AI1 The spin part of 

O Knight shift for the field H = 0.65T parallel to ab plane does not change down to 15mK, much below 
T C (H) « 1.2AT13. Since the experiment has been performed on S^RuOi in a clean limit this constitues 
the evidence for the triplet spin pairing with spin of paires lying in ab plane. On the other hand the 
measurements of the basal plane upper critical field down to 0.2J*T shows no tendency to divergency of 
H c2 {T^A ( see also results of measurements on less perfect crystalsEJ). Moreover, the search for reentrance 
of superconductivity under magnetic fields up to 33T and temperatures down to 50 mK has given the 
negative resultsE3. The basal plane upper critical field saturates at low temperatures at 1.5T which is 
well below the paramagnetic limit H p ss 2.8T and roughly corresponds to quasiclassical upper critical 
field value 

H c2 = ° , (1) 

where $o is the flux quantum and £ a f,, £ c are the coherence lengths in basal plane and along the c axis 
correspondingly. P 

These observations being certainly in contradiction with theoryfcl stimulate us to reinvestigate the 
problem of upper critical field in quasi-two-dimensional superconductors taking into account scattering 
of quasiparticles on the impurities. The main goal is to find a limits for crystal purity when one can hope 
to see the low temperature upper critical field divergency. In a few words the result can be described as 
follows. In a pure crystal the H C 2 divergency starts to be developed when the thermal coherence length 
£(T) = vp/2wT begins to be larger than the "cyclotron radius" of the quasiparticles orbits R C (H) = 
Vp/ui c . Here Vp is the basal plane Fermi velocity, to c — eHvpd/c is the "cyclotron frequency", d is the 
distance between conducting layers. The impurities does not prevent the upper critical field divergency 
if at the temperature determined by equality 

£(T) = R c {H c2 ) (2) 

the quasiparticles mean free path I is still larger than £ (T) . Otherwise the upper critical field is saturated 
at temperatures below 

and its value is roughly described by formula ([!]). 
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The paper has following structure. The general formalism for the upper critical field problem in layered 
conventional and unconventional superconductors is introduced in the next Section. By the way of the 
derivation the several simplifications have been used: (i) The equations are written for axially symmetric 
crystal where the only one-dimentional, even or odd in respect to the reflections in (H, crystal axis) 
plane superconducting states are realized; (ii) Among them the only even superconducting states are 
considered; (iii) The equations are derived with but solved in neglect Pauli paramagnetic interaction; 
(iv)Thc equations are solved for only unconventional superconducting states with zero anomalous Green 
function self-energy. The analytical solution of the equations accompanied by the discussion of limits 
of crystal purity sufficient for the upper critical field saturation at low temperature is containned in the 
third Section. 



II. THE ORDER PARAMETER EQUATIONS 



The electron spectrum of a layered crytals obeys basal plain anizotropy, z-dependent corrugation of 
the Fermi surface and several bands in general case. It seems however unimportant for our purposes to 
take into account all these complifications. So, we shall consider a metal with electron spectrum 

e(p) = +Pl) - 2tcos Pz d, 

t«e F = ^, —<P,<- d (4) 

in the magnetic field H = (0, H, 0), A = (0,0,— if a;) parallel to the conducting layers with distance d 
between them. In absence of impurities the normal state electron Green function G^^^iPyiPz, %, x') is 
obtained as the result of solution of the equation 



1 / d 2 2 \ . „. / , u> c x 



2m \ dx 2 



p + 2t cos p z d + OfJLeH + jl 



V F 



G U!(T (py,p z ,x,x')=6(x-x'), (5) 



where uj n — irT(2n + 1) is the Matsubara frequency, uj c = evpdH/c, vf is the Fermi velocity, h = 1, 
fi is the chemical potential, /i e is a magnetic moment of an electron in a crystal. To use the diagonal 
shape of the Green function matrix we have chosen the y direction as the spin quantization axis, such 
that a = ±1. 

In a layered crystal with a singlet Cooper pairing the superconducting states obey the following order 
parameters 

A s = iA s (p Xl p y ,p z ,x)a y , (6) 

Here, a y is the Pauli matrix. For a triplet superconductivity we shall limit ourselves by the consideration 
of the so called equal spin pairing states with spin lying in the plane of the conducting layers. In neglect 
of small effects of spontaneous magnetismtfl a vector wave function of such the states is 

d = zA t (p x ,p y ,p z ,x). (7) 

As we have put the spin quantization axis along y direction we will use the corresponding basis of Pauli 
matrices if = (<j y . & z ,<j x ). So the order parameter for triplet pairing state in our case is 

A* = i(do>)ay = -A t (p x ,p y ,p z ,x)a z . (8) 

The order parameter function is representedEl as a linear combination of the basis functions vpi(cf>,p z ) 
of one of the irreducible representations of crystal symmetry group 

A^{p Xl p y ,p z ,x) = ^, Pz )A s -\{x). (9) 

Here, <j) is the angle between basal plane vector of momentum p|| and magnetic field H || y. There are 
only one and two dimensional representations (i = 1, ...d; d = 1, 2) in the crystals with axially symmetric 
spectrum (Q). 
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The crystal with axial symmetry under magnetic field lying in the basal plane obeys the symmetry in 
respect of reflections in plane where the vectors of magnetic field and crystal axis lie, that is (y, z) plane 
in our case 1 . The two components of vector basis functions (ipi(4>,p z ),ijj2{(f>,Pz)) can be always chosen 
such that each of them will have definite (even or odd) and at the same time mutually opposite parity 
in respect to reflections in (y, z) plane. As the consequence the set of the order parameter equations for 
two component superconductivity splits on two independent equation for each component of the order 
parameter. One of them corresponds to the higher value of the upper critical field. Thus we always 
deal with one-component superconductivity with definite parity. The simplest examples of even functios 
ip(4>,p z ) are: 1, sin 2 <fi — cos 2 <fi for singlet pairing and sm(p z d), cos 4> for triplet pairing. As an examples 
of odd states one can pointed out on sin(p^c?) sin <fi for singlet pairing and sin <p for the triplet pairing. 

The upper critical field is found from the equation on the order parameter which have the different 
shape for even and odd superconducting states. For determiness we shall consider just the even order 
parameter states. In this case the order parameter equation for a clean superconductor with singlet 
pairing is 

Tl/d 

A s (<f >) p z ,x)=gJ dx' j^jL J d A^ VzW ^^ z) 

-Tr/d 

X T Y.\ E G^Av^v' z ^^)G-^-o(-v' y -v'^^)^^,v' z ^ (10) 

n a— ±1 

and for the triplet pairing 

Tt/d 

A t ( ( f,,p z ,x)=gjdx' J ^j(<b, Pz )rW,p' z ) 

-n/d 

x T Y,\ E G^Av'y^ z ^^)G^A-v' v ,-v' z ^^)^{4>',v' z ,x') (ii) 

n a— ±1 

In presence of the impurities the order parameter of a superconducting singlet pairing state 
A s = iA s (4>,p z ,x)a y acquires a self energy part 

S s K,i) = iA s Un (x)a y + A~l n (x)a x (12) 

consisting of singlet and triplet components ^| 

For a triplet pairing states with an order parameter A* = — A t ((j),p z , x)a z the corresponding self energy 
part 

t\u n ,x) = -Al n {x)a z + iAl n {x)a (13) 

consists of two different triplet components. Here &q is the two dimensional unit matrix. The singlet 
component of self energy part is absent for chosen basal plane orientation of magnetic field and equal 
spin triplet pairing with spin directions parallel to conducting layers. r— , 
The self-consistency equations for the order parameters and self energy parts have the form (seet3) 

-Tr/d 



1 For the uniaxial crystals with hexagonal or tetragonal symmetry this property takes place only for particular 
directions of magnetic field-ki the basal plane. 
2 Compare with the papeitj, where a similar theory in frame of quasiclasical approach has been developed. 
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n 



x GZ«(p' y ,p' z , x, x')G% n (-p' y , -p' z ,x, x')[A^ a(3 ^',p' zl x') + E*-* oi8 (w n ,a/)], (14) 

where V£ q ,a/'* = 5 3 s '' Pag 3 ' 1 + \^(<f>,p z )i>*(<l)',p' z )/2, .9 s = iS^ <?* = ™<7 Z ; 
G^ = G 0n ,i(ffo X7 + ^ A7 )/2 + Ga B ,_i(ff A7 - ^/ 7 )/2. 

Here tD n = o; n + Tsignuj n , T = mnu 2 fid and it, n are impurity potential amplitude and concentration. 
We will use also a quasiparticle life time r and a mean free path I introduced by means T = l/2r = vf/21. 

The equations (^) ace. obtained in frame of procedure of averaging over an impurity positions. As it 
was shown in the paperllll one can use the field independent value of T so long 

- > (15) 



The system of "the scalar" self-consistency equations for the singlet pairing states following of equations 
(0) is 



A s ^,p z ,x)=gfdx' J^L J d A^ VzW ^^ z) 

-7T /d 



r Ej E G^, CT ( P ;,^,:r,:rOG_^,- CT (-p^ (16) 



x ^ , 

2 

n cr— ±1 



Txjd 



A Q n ( x )= nu2 J dx ' J J E G ^ n AVv^Vz^,x')G^ n ^ a {-p vl -p zl x 1 x') 

-■Kid a=±1 

x [A s (^y) + A?>') + ^>')L (17) 



7r/d 

A L ( x ) = nu2 J dx ' J "7^ y 7^ E G Zn,<r(Pv>Pz> x > x ') G -u,„,-*(-Py,-Pz,X,x') 

-Tr/d CT=±1 

x [a(A s (^p z ,x') + A?J^)) + A^(a:')]- (18) 
For the triplet pairing case the corresponding equations are 

Tl/d 

, f d Py f dp' z 



A\<j>,p z ,x)=g J dx' J J |^(<^)W,P*) 

-7r/d 

xT E^ E G^^(^,p;,a ; ,xOG^„, ff (-p;,-^,x,xO[A t (^^p;,a ; ') + A^(a ; ')-^aA^(x')], (19) 
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7r/d 



K -*ld ^ CT=±1 

x {A^p^x^+AUx^-iaAUxt)}, (20) 
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ir/d 

Ai„(x)=nu 2 I dx' I -P- I ~!r~\ G ^^APy,Pz,x,x')G- Gjn ^(~p y ,~p z ,x,x') 

u ^ <r=±l 
—rr/d 



x (ia(A\cj>,p z ,x') + AUx')) + AUx')]. 



(21) 



Taking into account that in common the paramagnetic limit of superconductivity H p is much higher 
than an orbital upper critical field we shall rest the general problem of influence of paramagnetism 
on superconductivity for a future investigations. In the absence of the Pauli paramagnetic interaction 
the equations (|T^)-( 21 ) are greatly simplified and we have the system of two equations with equivalent 
structure for singlet and triplet pairing states 



Tr/d 



A" 



p zi x)=g I dx> I %L f ^fapjrWJ,) 



2tt 



2tt 



-Tr/d 



x Tj2G& n (p' y ,p' z ,x, x')G- Qn (-p' y , -p' z ,x, x')(A^(^,p' z , x 1 ) + Ajf(x')), 



(22) 



-Tr/d 

x G Qn (p y ,p z , x, x')G^ Qn (-p y , -p„x, x'){A^{^p Zl x') + A s / n (x')). (23) 

Below we shall omit the supercripts s,t using the common notation A s, *(0,p z , x) = A(<f>,p z ,x), 

A s ~ (x) — A^ n (x) for the order parameters and the self energy parts both in singlet and in triplet case. 

On this stage it is useful to note that the normal metal electron Green function which we should find 
as a solution of the equation (|^) depends of p y only through its square. Hence for some unconventional 
superconducting phases like A((j),p z ,x) — v^sin 2 cf) — cos 2 <fi)A(x) for singlet pairing or A((j),p z ,x) — 
V / 2cos0A(cc) for triplet pairing, the following property takes place 

We will discuss further the only unconventional superconducting states obeying the property (p4|). For 
such the states the self energy part is equal to zero and we deal only with the order parameter equation. 

Tr/d 

A(d>,p z ,x)= g J dx' J ^ J ^{^ PzW {4>\p' z ) 

-rr/d 

x T^G^p^p^x^^G^i-p^-p^x^^A^^p^x'). (25) 

n 

The equality ( p4| ) is not valid for conventional superconducting state as well for many unconventional 
superconducting states where, as the consequence, there are nonzero selfenergy parts. In isotropic conven- 
tional superconducting state it prevents a supression of the superconductivity by the ordinary impurities. 
For unconventional superconducting states the presence of the self energy leads just to the mathematical 
complifications and does not change qualitatively the main results. 

The only difference of ( pB] ) from the pure case consists of change oi„ — > uj n and one can use the expresion 
for the normal metal electron Green function found in the paperEJ. It has nonzero value in the regions 
determined by the inequality u> n (x — x±) > for positive value of x component of electron momentum 
(a = 1) and by the inequality to n {x — x±) < for negative values of x component momentum ( a = — 1), 
where it is defined as 
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—isignujn 
Dp sin < 



x exp 



Here A = 4t/u) c . 

This expression is valid under the condition 



■ exp 



au) n (x - xi) 



f iaX . 


LO c (x - X\) 






cos 


^ sm 


2v F 





Vzd- 



vf sm0 
to c (x + X\) 
2^ 



e-np[iapF (x — xi) sin < 



(26) 



I sin 0| > \l — . 



(27) 



The disregard of the small intervals of the angle (f> where | sin0| < (t/ep) 1 ^ 2 means that the only open 
trajectories of quasiparticles on the Fermi surface are taken into account f\ They give the main singular 
part to the kernels of the integral equation (P5|) . 

The substitution of (|26|) to the equation (|27|) gives after summation over the Matsubara frequency 



A(x)= g 



\x — xi\>a sin4 

2iX 

x exp < — — - sm 



dx\ 



2itvf sine 



d(p z d) 



2tt 



tp*((j),p z )ip{cj>,p z )- 



2-kT exp 


\x—x\ 


I sin cp 


sinh 


"27rT|a;-a;i|" 
vf sin <p 





uj c (x - Xl) 

2^p 



sin(p z <i) sin 



LU c (x + Xl) 



2vi 



(28) 



Here g = am/And is dimensionless constant of pairing interaction, a is the small distances cutoff. The 
equation ( ps| ) is the basic equation of the paper. The properties of its solution we shall discuss in the 
next Section. 



III. THE UPPER CRITICAL FIELD 

Let us make an unessential simplification and consider p z independent superconducting states deter- 
mined by functions normalized as follows 



In this case one can perform the integration over p z in (E 



A(x)=g 



dxi 



2nvp sin 

x — a;i|>a sm0 

uj c (x - Xl) 



r ^*(0>(0> 



2wT exp 


\x—x\ | 


I sin (j) 


sinh 


' 2-kT\x-x x \~ 
vf sin <f> 





X 2q < — r sm 



2v F 



sm 



uj c (x + Xi) 

2m? 



A(xi 



(29) 



Here 2o(...) is the Bessel function. For the further purposes it is not important to work with the equation 
( p8| ) or with (^9|) and for the determiness we will operate with the latter. 
The equation (p9|) in the absence of a magnetic field 



3 The numerical calculation taking into account both open-,and slosed trajectories, performed for clean case in 
the papenU, just confirms the qualitative results of the articlea have been obtained in neglect of closed trajectories. 
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oo 

2naT 

determines of the critical temperature T c , which is expressed from here through the critical temperature 
T c q in a perfect crystal without impurities I = oo 

T c0 = ^exp(-l) (31) 



na \ g 
by means of well known relation 

h £-*(i)-*(i + ss5;)- < 32 > 

For the critical temperatures T c ~ T c q the suppresion of the superconductivity by the impurities is: 

T c = T c0 -^- (33) 

One can also point out the condition of complete suppresion of the superconductivity 
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(34) 



where In 7 — C = 0, 577... is the Euler constant. 

The behavior of the upper critical field is determined by the relationship of the three spacial scales: 
vf/2ttT, I and vp/uj c . For the temperatures near to the zero field critical temperature T ~ T C (H = 0) 
the upper critical field H C 2(T), or us C 2(T) = eH C 2(T)vpd/c tends to zero and the inequality 

min {^f i } < ^5ry (35) 

always presents . 

At the temperature decrease the upper critical field increase, but so long the inequality (^35j) takes place 
the essential interval of integration over [x — x%) in ( p9|) determined by the min(vp/2iTT, I). Hence, one 
can use the smallness of (x — xi)vf/oj c in the kernel of the equation (p9j): 



Ah) = ,/ / ,/,, / . r (.>„■(,„ ,U(X ^ xJ ^Uc*. A()|) 



I'kvf sin 



|x — a;i|>a sin0 



27rT exp 




/ sin 


sinh 


"27rT|a:— a?i|~ 
vjr sin 





lysine) 



(36) 



The solution of this equation gives the correct value of w C 2(T) so long the inequality (|35J) is truth. By the 
substitution of the uj C 2 (T) obtained from this equation to the inequality ([55|) one can establish the upper 
limit of the superconductor purity at which the H c2 (T) found from this equation represent the correct 
value of the upper critical field up to zero temterature. 

For the pure enough samples and low enough temperatures it can be turn out that the opposite to the 
( |35| ) relationship 

vp . ( Vf 

< mm 



W c2 (T) 



4 Near T c there is also formal solution of ( p9| ) with u> C 2 > t ( seeS). This solution, however, is related to the region 
of magnetic fields, out the region of validity of the present theory (|15[). 
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breaking the correctness of transfering from ( |2i| ) to fl36[ ) is realized. It should be noted that ultraclean 
case demands a special investigation (see condition ([15]) and discussion in the paperEZI). One can claim 
however, that at the temperatures 

Vf t u; c (T) 
2tt2 2tt 

the magnetic field dependence of the critical temperature in the equation ((2^) starts disappear or, in 
other words, the tendency to the divergency of the upper critical field pointed out in the papercl appears. 

To solve of the equation (|3^) at arbitrary temperature and purity is possible only numerically. Here 
we shall discuss the case when it allows the analitical solution. If the legth scale £, on which the function 
A(x) is lokalized, is larger than the essential distance of integration over (x — x\) 

i > min {^ l }> ( 38 ) 

then one can expand A (x±) w A(x) + A'(x)(x— x%) + A"(x)(x — xi) 2 /2 under the integral in (|36|). Taking 
into consideration that under this condition the argument of Bessel function turns to be small even on 
the upper boundary of effective interval of the integration over (x — x\): 

tu c2 (T)£ min{^ r: 1} _ t min{j^, 1} ^ ^ 
one can also expand Bessel function Tq(x) sa 1 — x 2 /4. As the result we get the differential equation 

2 



^-K^) + ^))*'>--^(&)"^> + '<">( 



tu> c2 (T)x 



ttvfT 



A(x), 
(40) 

where a = (2"kTt)~ 1 , 

oo 

1(a) = / — -exp(-az)^4^ (41) 
J smhz (2n + 1 + a) 6 

n n=0 v ' 



G+ = / — ^(^^(^sin^. (42) 



7T 



(I 



In pure case a w a c = (2itT c t) 1 <C 1 the inequality ( |38| ) is valid only in vicinity of the critical 
temperature. Puting T — T c in the right hand side and taking its lowest eigen value we obtain 

, T c0 (1 1 \ /1\ v ^/(a e )^ e2 (T) 



Summing this equation with equation ( |32j ) and leaving only linear on T — T c and on the impurity con- 
centration terms we get 

Wc2 (T) = ^H c2 (T) = -^=- (T C0 - /3f ) (T c - T). (44) 



Here the coefficient 
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shows that, the slop of H C 2(T) at T — T c decreases with the increase of the impurity concentration 
somewhat slower than T c itself (see eqn. (^)) . 

The equation ( |36| ) does not contain any divergency of H C 2 (T) and the expression ( [44| ) which is valid 
in Ginzburg-Landau region can be used at arbitrary temperature as the estimate of upper critical field 
from above. Hence, to establish the limits of a sample purity, at which the equation (36) works, one may 
substitute ( p4|) at zero temperature into the inequality (fiEl). Omitting the numerical factor of the order 
of unity we have 

I < JfrU, (46) 

where £ a b = vf/2ttT c is the basal plane coherence length. We see, that there is the good reserve in sample 
purity in the limits of which one can not expect low temperature divergency of the upper critical field. In 
the Sr2RuOi the mean free path should be approximately 10 times larger than the basal plane coherence 
length £ to go out limit fl46|). 

Let us consider now the dirty case: T c <C T c o, a>l and 1(a) sa a~ 2 allowin g an alytical solution for 
H C 2(T) at arbitrary temperature. Taking the lowest eigen value of the equation (|40|) we get 

Too , A, 1 \ , /l 



111 ~F " ^ I 2 + 4^T ) + ^ [ 2 ) = V2G^^ca(T). (47) 
Summing of this equation with equation ( |32| ) yields 

"« {T) = (48) 

This expression is correct at any temperature. One can rewrite it approximatively in more simple form 

Wc2 (T) = -^L(T c 2 -T 2 ), (49) 



where 



T c = ^n^Y /2 . (50) 



IV. CONCLUSION 



The system of linear integral equations for the order parameter of conventional and unconventional 
superconducting state in a layered crystals with impurities under magnetic field parallel to the conducting 
layers is derived. It is shown that so long the purity of the samples does not exceed high enough level 
( [46| ) there is no tendency to the low temperature divergency of the upper critical field. The analytical 
solution of the equations in the clean (Ginzburg-Landau region) and dirty (arbitrary temperature) limits 
are presented. 
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